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Abstract 



Recently, Anderson et al. [Ij published an empirical prediction formula for the so far unex- 

I plained parts of the velocity changes of spacecrafts during Earth flybys. In the framework of a 

■ perturbational approach, we show that there is no velocity-independent force field of the Earth 

, - in addition to its Newtonian gravity field - that is to reproduce this formula. However, we 

I give examples for fields modeling exactly the flyby anomaly which are quadratic functions of 

^ ' the velocity of the spacecraft. 

CN ■ 1 Introduction 



CJ I The flyby anomaly is an unexplained frequency change in the radio Doppler data during an 

Earth flyby of a spacecraft. Interpreted as a "real" effect, this corresponds to an anomalous ve- 
locity change by a few mm/s as we can see, for example, in Table 2 in the paper of Lammerzahl, 
Preuss and Dittus [2|, which gives an overview of this phenomenon. 

Recently, Anderson et al. |T] published an empirical prediction formula (eq. (2) in [IJ, eq. 
^ ■ (1) in section 2 of the present paper), which accurately reproduces the measured anomalous 

^\ I velocity changes for all of the up to now six flybys for which sufficient data are available. 

Despite the fact that this formula is not related to any explanation of the phenomenon, it is a 
valuable progress for two reasons. Firstly, it makes sure that the anomalous velocity changes 
, are not only significant, but there is a clear regularity in the data that rules out any exclusively 

I statistically disturbing influences on the spacecraft orbits. Secondly, the special form of the 

discovered regularity enables us to systematically draw conclusions about the form of the 
unknown forces. Especially, if we assume that there is some unknown force field generated by 
the Earth (in addition to its Newtonian gravity field) and acting on the spacecrafts, we are no 
^ ' longer in the situation to guess the form of this field and then check whether it explains the 

measured anomalous velocity changes in all six cases. Instead, we can try to find force fields 
reproducing this formula by much more efficient mathematical methods. Certainly, this has to 
be done under the assumption that the prediction formula applies to all possible Earth flybys 
(at least in a wide range of the relevant orbital parameters) . Finally, we can draw conclusions 
about the general properties of all such flelds. 

In the following section the latter will be done, resulting in strong restrictions for those 
fields. In particular, we prove that any force which depends only on the spacecraft's position in 
space and not on the spacecraft's velocity (like the Newtonian gravity) is strictly ruled out. On 
the other hand, section 3 demonstrates there is a wide variety of force fields being quadratic 
functions of the spacecraft's velocity and modelling exactly the fiyby anomaly. In section 4 we 
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discuss our results in view of the consequences for the possible nature of the unknown kind of 
interaction and for the energy transfer process. 



2 Exclusion of velocity-independent force fields 

The approximation formula for the anomalous part AVoo of the velocity change given by 
Anderson et al. [l] reads 

—— = -— = K{cos 6i - cos 60) , (1) 

where Voo denotes the hyperbolic excess velocity of the spacecraft; E and A£^, the total specific 
energy and the anomalous part of its change respectively, and bi and the declinations of 
the incoming and outgoing asymptotic velocity vectors. (The declinations, as well as Voo and 
are defined as osculating Kepler orbital elements.) IT is a dimensionless constant which is 
in good agreement with twice the special relativistic /5-factor of the Earth's spin rotation at 
the equator: 

K = 2i^^^ « 3.099 • 10-'^ (2) 
c 

{tOE^ Re'- angular rotational velocity and mean radius of the Earth. Even if we regard values 
specific for the Earth, hereafter we will not subscript the Earth's mass and thus it is Me = M.) 

Our goal is to answer the following question: What is the mathematical form of a (classi- 
cal,) stationary force field acting on a spacecraft (additionally to the Newtonian gravity field) 
that reproduces formula ([T])? Since the anomalous accelerations necessary to yield the mea- 
sured anomalous velocity changes are only of the order 10~^ of the Newtonian gravitational 
acceleration [2113]) it seems to be sufficient to treat the additional forces as small perturbations 
along the classical Kepler trajectories of the spacecrafts. In this section, we will restrict our 
considerations to forces that only depend on the spacecraft's position in space but not on the 
spacecraft's velocity. According to ([T]), such a field F must yield a total specific energy change 

00 

A^= y" F {r {t)) ■ V {t) dt = 2KE {cos 5i- cos 5 o) (3) 

—00 

for any hyperbolian Kepler orbit. The integration is to be taken along the unperturbed tra- 
jectory r{t) with velocity vector V = r. (Here and in the following E, AE, F etc. denote 
the specific quantities, what means that, as in [1], these values are divided by the spacecraft's 
mass. Therefore, i. e., in accordance with the structure of ([T]) and the force has been 
replaced by the acceleration.) 

The most general expression for the field F = F{r, a, 5) with r = \r\ is 

F{r,6,a) = fi{r,6,a)er + f2{r,6,a)es + f3{r,6,a)ea , (4) 

where e^, es and Sa are the unit vectors belonging to the radial and the usual equatorial 
geocentric coordinate directions respectively (they form an orthonormal basis at each point of 
the three-dimensional space, and, of course, a is the right ascension.) In view of ([T]) and 
the source of the additional field should be in a today unknown relation to the rotation of the 
Earth. Concerning the question whether there is a solution to equation ^ at all, we will show 
now that without loss of generality we may assume that /2 = /s = and fi{r,5,a) = 
with some symmetric function /, appealing to the symmetry properties of formula ([T]). (Note, 
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however, that to any sokition F there are more solutions without the symmetry properties 
of ([T]). One possibihty to find such solutions is to add the gradient of a potential which does 
not share the symmetries and converges to a constant direction-independent value in the limit 
r — > oo.) 

Step 1: 

Due to the rotational symmetry of ([T]), with respect to the Earth's axis, for any solution F 
there are new solutions, obtained by shifting the a-coordinate by an arbitrary constant value. 
By means of superposition of "shifted" solutions of the equation ([3]) (which is linear in F), we 
can completely eliminate the a-dependence of the functions /i,/2 and /s (within the limit of 
the superposition of an infinite number of solutions) . 

Step 2: 

Because of the invariance of ([T]) under the transformation a — > — a (transforming ■ V to 
—Ba • V), we know that for any solution F{r, 6) = fi{r, 6)er + /2(?", S)es + fs^r, 6)ea the field 
/i(r, 6)er + f2{r, 6)es — f3{r, 6)ea is also a solution. Half the sum of both solutions has no Cq, 
component; so without loss of generality we may henceforth set /a = 0. 

Step 3: 

Similarly, the mirror symmetry of ([T]) with respect to the Earth's equatorial plane implies that 
with any solution F{r,5) = fi{r,6)er + f2{r,S)es the field fi{r, —6)er — f2{r, —6)es is also a 
solution. Half the sum of both solutions has a symmetrical part and an antisymmetrical 
es part; so without loss of generality we may assume that fi{r,6) is symmetric in 6 whereas 
/2(r, (5) is antisymmetric in 6. 

Step 4: 

Formula ([T]) reveals another symmetry property, namely a scaling symmetry: Let us be given 
an arbitrary hyperbolic orbit parametrized by the polar angle (/? of the polar coordinates in 
the orbital plane, i. e., r{ip) with modulus 

r{<f) = ^ , (5) 

1 + e cos ^ ^ 

where p and e are the semi-latus rectum and the numerical exentricity, respectively. If we 
enlarge this orbit by a "blow-up" factor b G (0, oo), we get a new orbit br{ip) with total energy 
y instead of E. The other terms on the right-hand side of equation Q remain unchanged 
under this transformation; so the anomalous energy change AS transforms also with the factor 
^, which yields the scaling property 

b j F{br) ■dr = ^j F{r) ■ dr (6) 

with the factor b on the left-hand side coming from the streching of the integration path (both 
integrals are understood to go along the unenlarged (6 = 1) path in space). Since this scaling 
property holds true for all hyperbolic Kepler orbits, the effective part of the field F contains 
its r dependence only in form of a factor whereas other powers of r do not contribute to 
the net energy change AE. (Formally, one can work this out by a Laurent development of F 
with respect to r, followed by a comparison of the powers of b.) 

Consequently, if there is any solution of equation (jSj) , we can get a new solution by discard- 
ing all parts but the ^ part (irrespective whether the original solution is of the "long" form 
(m or of a reduced form after the steps 1, 2 or 3). Hence, going on from step 3, without loss 
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of generality we assume that F is of the form 

where /i is a symmetric and /2 an antisymmetric function. 
Step 5: 

Finally, we may even omit the eg part of ([7]), as can be seen by the following "radialization" 
procedure: Replace F{r, 6) by F{r, 6) + S/^^, with $ being an antiderivative of —f2- (Because 
of lim = 0, the additional term does not contribute to AE.) In our spherical coordinates, 

r— >oo ^ 

the gradient term reads = — ^ffie^ + es and further, together with ([7]) and = — /2, 

it gives F{r, 5) + = ^^^^^^^^^^ e^. Thus, if there is any solution of the form ([7]), there will 

be a solution of the form 

Fir,S) = ^er, (8) 

too (with / = /i — $). Since /2 is assumed to be antisymmetric, <i> is a symmetric function, 
which guarantees that the symmetry of the part in 6 is preserved under this manipulation. 
(Alternatively, by an analogous procedure, one can transform the part of ([7]) to zero, leaving 
a purely tangential field.) 

Now, to prove the velocity independence of any force field wrong, we assume there is a 
velocity-independent force field of the form ([8j) that reproduces equation ([3]) for any hyperbolic 
Kepler orbit and lead the assumption to a contradiction. It turns out, that it is sufficient to 
deal with a specific class of polar orbits, namely orbits which go from south to north in such 
a way that on the incoming branch and during the passage of the perigee 6 is positive and on 
the outgoing branch the spacecraft passes the pole (5 = ^). In the following, we will restrict 
our considerations to these polar flybys, for those there is 

COS Oi — COS do = —I sm ( — — ) sm ( — — ) = I- cos Op (9j 



e 

with bp being the geocentric declination for the closest-approach location. By that we have 
applied the general relation - = sin ^ between the numerical eccentricity e and the deflection 
angle A of a hyperbola. Remarkably, we see that ujeRe cos dp = Vp, where Vp defines the 
tangential rotation velocity of the Earth at the perigee. Hence, the specific energy change 
depends directly on this velocity at the perigee, what will be discussed later on more explicitly 
in Section 4. 

According to £^ = = ^{e — l)^r^ where is the value of the radius coordinate r at 
the perigee, equations ([3]) and Q lead to 

AE = 2K^^^ cos 6p . (10) 

£ Vp 

On the other hand, we can express the energy change by inserting a field of the form ([8]) into 
equation ([3]) which gives 

oo 

AE = KGM I ^rdt. (11) 
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Introducing the polar angle ip for the parametrization of the hyperbola (see ([5])), we can write 



f d /1\ d (l + ecosip\ el ^ . 

sin(/9(/9. (12) 



dt \r J dt \ {1 + e)rp J l + evp 
With that the comparision of the right-hand sides of equations (jlOp and (llip result in 



£2 



2cos(5p = j f {5) sin dip . (13) 



For e — > oo, whereby the hyperbolas degenerate to straight lines, we get 

ZL 
2 

2cos(5p= j /((5) sin d(/3 . (14) 

TT 

2 

For our special polar flybys this limit yields ip = 5 — 5p ioi — '^ < ^ < ^ ~ ^'^d ip = ir — 6 — Sp 
for — 5p < ip < ^. Consequently, we should divide the integration according to the borders 

2 cos 6p= J /(5)sin(<5-5p)d<5- j f (6) sinin - 5 - 6p) d5 (15) 



cos (5p 



f 2 2 "PN 2 

11 / /Wcos5d(5. (16) 

\-f+<5p f / -f+^p 

Since / can assumed to be a symmetric function (see Step 5 above) the part 6 £ [~f + ^pj f ~ ^p] 
of the first integral vanishes. After division by 2cos5p, equations (fT5|l and (fT6]) reduce to 

f f-5p 
1= j f{5)sm6d6 - tan 6p j f{6)cos6d6. (17) 



Furthermore, for (5p — > the right-hand side of the equation vanishes as well and the whole 
equation results in the inconsistency 1 = 0. Hence, we see that the assumption there is a 
velocity- independent force field cannot be correct. 

It may be noted that, because of the antisymmetry of AE as well as J V dt under the 
transformation t — > — t, terms with velocity dependence of odd-numbered power cannot con- 
tribute to the anomalous energy change AE. For scaling reasons the radius has the power of 
p — 4 when the velocity has the power of p. For any p > 2 the force field would not vanish for 
increasing radius. Therefore, other values for p than and 2 are out of the question. 



3 Velocity-dependent solutions 

Since we have seen that there is no velocity-independent force field reproducing equation (1) we 
will now put up with forces depending on the spacecraft 's velocity. We will see that intuitional 
transformations of ([1]) lead to a wide varity of fields which exactly model the flyby anomaly. 

An obvious idea to find velocity-dependent fields F{r, V) solving ([3|) is to regard the right- 
hand side as the difference of the values of an antiderivate of a suitable function, evaluated 
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at the "in" and "out" values of the integration variable, and "read off" the field from the 
integrand. In the following we can resort to the fact that the specific energy is a constant 
value along any undisturbed trajectory, but in order to get a field which depends only on the 
variables r and V (and in order to avoid obtaining a family of fields, parametrized by E), we 
have to express E hy r and V through E = ^V"^ — 
Formally, this approach gives for the present 

oo 
— oo 

but there is an ambiguity concerning the meaning of the 6 coordinate. Only their special 
values 6i and 60 are well defined by the directions of the asymptotic velocity vectors. However, 
since the asymptotic directions of the velocity vectors coincide with the asymptotic geocentric 
positions of the spacecrafts, we have the choice between the "velocity-space interpretation" 
and the "configuration-space interpretation" . The former interpretation leads via 



d . d I 1/2 



where is Vz := Vsm6, to quite "unnatural" fields with a complicated V dependence and 
components of V in the denominator, whereas the latter interpretation deals with 6 as a 
coordinate (a function) in configuration space and gives 

00 

AE = -K j {v^ - (V cos 5)-Vdt (20) 



-00 

00 



I ( -?-^^)sm<5e5- Vdt. (21) 



r 



-00 



The comparision with ^ shows us that 



(T/2 '2C'M\ 
--^jsinJe, (22) 

is a solution of Therefore, it reproduces exactly the formula ([T]) in the framework of our 
perturbational approach. 

Let us point out a remarkable property of this field. Considering the idea of interpreting ([1]) 
as an evaluation of an integral, it can easily be understood that the flyby anomaly affects only 
unbounded orbits. Since for each Kepler orbit the energy and thus V'^ — is a constant, a 
spacecraft is effectively subject to a field proportional to Vcos5 (cf. (f20]) ). which is therefore 
a conservative field (in the sense that the corresponding 1-form is not only closed, but also 
exact). Consequently, the integral ()20p vanishes for any bounded (elliptical or circular) Kepler 
orbit (in this case ()2ip comes from a loop integral, and t, of course, extends only over a 
finite interval). Concerning the Pioneer and flyby anomalies, this is in accordance with a 
conjecture often pronounced by Claus Lammerzahl (see, e.g., [2|, p. 97), namely that escape 
orbits behave qualitatively different from bounded ones, whereby the latter show no anomalous 
effects. (But note that, of course, our field (I22p gives a vanishing energy change in general only 
after integration over a complete cycle.) 

As easily can be seen, the force field (f22]l is discontinuous at 6 = Especially the 

divergence of the velocity-independent part goes to infinity at the axis of Earth. In connection 
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with this undesirable properties let us consider the term = — ^p-e,, — ^p-e^. Since this 

term multiplied with an arbitrary constant gives no contribution to the integral we may 
add -2wKGMV^ to our solution l\22\\ . where the "weight" w is an arbitrary real number. 
The resulting one-parameter family of solutions reads 

= -K [2{w -l)-^cos5er - i — -2w-^]sm6es] (23) 



where is 



' V 2GM 
and Fq = K [ — sin 5 es -\ — cos 5 ) . (24) 



(Note that for bounded orbits the property of vanishing total energy change mentioned above 
holds true as well for any w.) The special case w = {) corresponds to the " radialization" 
procedure in Step 5 of Section 2 and gives a field {Fq) with a continuous velocity-independent 
part. The divergence of this part vanishes everywhere. 

However, as we have analyzed, the velocity-dependent part of the field ()22p cannot be 



transformed in a similar way. Independent of the " radialization" , the problem with a diverging 
total field energy in the limit r — oo of that part remains unsolved. 



4 Discussion 

We have analyzed the consequences of the assumptions that the approximation formula given 
by Anderson et al. p!] is valid for all possible spacecraft flyby manoeuvers on the Earth and 
that the flyby anomaly is a "real" effect which can be described by a force field belonging to 
some unknown interaction between the Earth and the spacecrafts. (The interesting question, 
whether for flybys on other celestial bodies the coefficient K in equation ([2]) is to be replaced 
by twice the (3 factor of the body's spin rotation velocity at its equator, is beyond the scope 
of the present paper.) By a "real" effect we mean that the conclusion from the radio Doppler 
data to the anomalous velocity changes is correct and that the anomaly is not (totally or to a 
considerable part) due to unexplained influences on the propagation of the radio signals or to 
any kind of mismodelling in the domain of conventional celestial mechanics or other applied 
branches of physics. Since the mass of the spacecraft does not enter equation ([T]), it seems that 
the velocity changes are caused by unknown gravity-like forces, which, in a classical approach, 
are to be described by a force field. 

As we have seen, such a force field inevitably contains velocity-dependent terms in order to 
be in accordance with ([T]). Equation ([T]), however, has been published by Anderson et al. [Ij as 
an "empirical prediction formula" which may have only approximative character. It remains 
the possibility that for extreme values of the relevant orbital parameters this formula loses its 
applicability. To us, it seems reasonable to assume that the symmetry properties we utilized 
in steps 1 to 3 of section 2 are of general validity. Concerning step 4, one may argue that 
there may appear deviations for very large perigee distances. As long as such deviations are 
restricted to the region beyond the observed parts of the asymptotic branches of the hyperbolic 
orbits, this does not effect our no-go result. Let us also note that the no-go result does not 
depend on the validity of the approximation formula ([1]) for arbitrarily large excentricities 
e: Coming from ()13p . one can complete the proof without performing the limit e — > oo, but 
compare the powers of e in (fT3]) . However, due to the e dependence of the borders of (fT3|) 
(which diappears in the limit e — > oo), this calculation turns out to be much more involved. 

In the proof of our no-go result in section 2 we learned that the total energy change 
/\E remains bounded in the limit of infinitely large excentricity e, see (|10|) . (Contrary to a 
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first impression, this is not inconsistent with the fact that, according to formula ([T]), AE is 
proportional to the energy E which diverges in the limit e — > oo, since in this limit the cosine- 
difference factor goes to zero.) However, despite this property there is the problem of the 
temporally resolved amount of energy transfer in the course of the passage, which has already 
been studied by Anderson, Campbell and Nieto [4]. In the case of the field (f22l) . the energy 
transfer integrated over the incoming branch of the hyperbola as well as the transfer on the 
outgoing branch diverges in the limit e ^ oo. Only their sum remains bounded. (Instead of 
calculating the integrals belonging to the two branches, we can consider the situation within 
the picture of an effective potential proportional to E cos 6, where E is to be regarded as a 
constant. In this picture, it becomes clear that in the limit e ^ oo there are infinitely high 
potential differences to stride through.) 

As a byproduct of the relation ([9]) we see that for the special class of orbits, considered 
in section 2, the right-hand side of formula ([T]) reduces to K [cos Si — cos 6o) = 4^^^^cos(5p. 
Hence, the relative velocity change is proportional to the velocity Vp of the Earth's surface 
under the perigee. We want to point out that this interesting relation, which may feign that the 
energy transfer is caused by the movement of the Earth 's surface directly below the spacecraft, 
cannot, as we have analyzed, be generalized to nonpolar orbits by fitting in a suitable factor 
taking into consideration the inclination. Even for polar orbits whose direction of the velocity 
vector swings over one of the poles {5 = ±§), this relation does not hold true. 

In section 3 we mentioned that there is, in principle, another possibility to "intrerpret" 
the 6 coordinate in the integral (jlSp . For the following reason it is expected that this way 
leads to fields with vastly different energy transfer properties: Whereas in our "configuration- 
space interpretation" 5 varies over more than half the celestial sphere, starting from 6i until 
it reaches its final value 5o, in the "velocity-space interpretation" the respective arclength is 
much shorter, namely given by the total deflection angle. Moreover, the " radialization" of the 
velocity-independent part of the field (end of Sect. 3) affects the energy transfer process. This 
may be relevant in view of an observation made by Anderson, Campbell and Nieto [^i The 
energy transfer process does not consist of monotonic increase (or decrease) of energy with 
time but is rather more complicated. 

Additionally to the energy transfer problem, which is due to the lack or very small number of 
data points between the asymptotic branches of the hyperbolas, all our efforts to infer the acting 
forces from the measured anomalous velocity changes suffer from the lack of any information 
about the force components perpendicular to the direction of motion of the spacecrafts. So we 
must realize that nowadays there is no unique solution of our problem (unless we make further 
strong assumptions), and we are restricted to more general statements, especially about the 
velocity-dependence. 

Let us note that, quite generally, any velocity-dependent force requires further forces due 
to the Galilei invariance (or Lorentz invariance). In our case, at any point in space we can 
transform the V'^ part of the field to zero by a Galilei transformation, whereby the whole 
force has to be invariant. That implies nontrivial transformation properties of the velocity- 
independent part of the field. This shows us that both parts may not be viewed as independent 
fields and that the velocity-independent part has a structure which differs substantially from 
the common Newtonian gravity field. 

In view of this close interlocking of both parts of the field, we should not overrate the 
divergence properties of the velocity-independent part analyzed in section 3. (After the "ra- 
dialization", this part can be interpreted as a gravity shielding which is most effective in the 
Earth's equatorial plane.) 
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